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ABSTRACT 

We derive the number density evolution of massive field galaxies in the redshift range 0.4 < z < 1.2 using 
the /T-band selected field galaxy sample from the Munich Near-IR Cluster Survey (MUNICS). We rely on spec- 
troscopically calibrated photometric redshifts to determine distances and absolute magnitudes in the rest-frame 
/T-band. To assign mass-to-light ratios, we use an approach which maximizes the stellar mass for any K-bwd lu- 
minosity at any redshift. We take the mass-to-light ratio, ./^ jL^, of a Simple Stellar Population (SSP) which is as 
old as the universe at the galaxy's redshift as a likely upper limit. This is the most extreme case of pure luminosity 
evolution and in a more realistic model ^/L^ will probably decrease faster with redshift due to increased star 
formation. We compute the number density of galaxies more massive than 2 x IO'^/j-^^q, 5 x lO'^/i-^^g, and 
1 X lO^^h^^^Q, finding that the integrated stellar mass function is roughly constant for the lowest mass limit and 
that it decreases with redshift by a factor of ^ 3 and by a factor of ^ 6 for the two higher mass limits, respectively. 
This finding is in qualitative agreement with models of hierarchical galaxy formation, which predict that the num- 
ber density of ~ M* objects is fairly constant while it decreases faster for more massive systems over the redshift 
range our data probe. 

Subject headings: surveys — cosmology: observations — galaxies: luminosity function, mass function — 
galaxies: evolution — galaxies: fundamental parameters 



1. INTRODUCTION 

Significant progress has been made in understanding galaxy 
evolution in recent years due to results of statistically signifi- 
cant and complete redshift surveys spanning the redshift range 
< z < 1 (Lilly et al. 1995; Ellis et al. 1996; Lin et al. 1997, 

1999) and due to detailed kinematic and stellar population stud- 
ies of small samples of galaxies (e.g. Kelson et al. 2001; van 
Dokkum et al. 2001; Vogt et al. 1996; Bender et al. 1996). On 
the theoretical side, semi-analytic models of galaxy formation 
(Kauffmann et al. 1993; Somerville & Primack 1999; Cole et al. 

2000) try to reproduce properties of galaxies in the local uni- 
verse and predict their evolution with redshift within the frame- 
work of a hierarchical CDM-dominated universe. 

The traditional observables used to characterize galaxies are 
unsuitable for studying the assembly history of galaxies, one of 
the most fundamental predictions of CDM models, since these 
observables may be transient. The best observable for this aim 
is, in principle, total mass, which is on the other hand very hard 
to measure. It has been argued that the best available surrogate 
accessible to direct observation is the near-IR /T-band luminos- 



ity of a galaxy which reflects the mass of the underlying stellar 
population and is least sensitive to bursts of star formation and 
dust extinction (Rix & Rieke 1993; Kauffmann & Chariot 1998; 
Brinchmann & Ellis 2000). The main uncertainty involved in 
the conversion of /T-band light to mass is due to the age of the 
population, amounting to only a factor of two in mass uncer- 
tainty for populations older than 3 Gyr. 

In this paper we use the /T-band selected field galaxy sample 
of the medium deep, wide-field Munich Near-IR Cluster Survey 
(MUNICS) to study the number density evolution of massive 
galaxies. See Drory et al. (2001) for details of the survey (MU- 
NICS 1 hereafter). Data presented here refer to ^ 5000 galaxies 
for which high quality VRIJK imaging covering 0.27 square 
degrees is available (see §2). The distances of the galaxies are 
derived from spectroscopically calibrated photometric redshifts 
(§3) and their masses from rest-frame /T-band luminosities (§4). 
We discuss the resulting integrated stellar mass functions at dif- 
ferent mass limits and their evolution with redshift (§5). 

We assume £2^ = 0.3, £2^ = 0.7 throughout this paper. We 
write Hubble's Constant as Hq ~ 100 h km s^' Mpc^', using 
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Fig. 1 . — Left panel: Comparison of photometric and spectroscopic redshifts for 310 objects in 5 survey patches (different symbols). Middle panel: The histogram 
of the redshift errors. The rms scatter is consistent with a Gaussian (dotted line: best-fit Gaussian) of a width a = 0.075 and an insignificant mean deviation from 
the unity relation of < A z >= —0.006. Right panel: The distribution of photometric redshifts (soUd histogram) and a best-fit analytic description (dotted fine) as 
well as the distribution of spectroscopic redshifts (dashed Une). 



h = 0.60 unless explicit dependence on h is given. 

2. THE GALAXY SAMPLE 

The galaxy sample used here is a subsample of the MUNICS 
survey, selected for best photometric homogeneity, good see- 
ing, and similar depth. Furthermore, in each of the remain- 
ing survey patches, areas close to the image borders in any 
passband, areas around bright stars, and regions suffering from 
blooming are excluded. The subsample covers 0.27 square de- 
grees in y (23.5), R (23.5), / (22.5), J (21.5), and K (19.5); the 
magnitudes are in the Vega system and refer to 50% complete- 
ness for point sources. 

Stars were identified following the procedure described in 
MUNICS 1 which relies on the combined information of PSF 
shape and colors in the J—K vs. V—I plane. The color cri- 
teria may also exclude z < 0.25 compact blue galaxies. Such 
galaxies are very unlikely to be present in the /T-selected sam- 
ple. Since we restrict our analysis to z > 0.4 this is anyway not 
a problem. 

The final catalog covers an area of 997.7 square arc minutes 
and contains 5132 galaxies. The fields included in this analysis 
are S2, S3f5-8, S5, S6, and S7f5-8. See Table 1 in MUNICS 1 
for nomenclature and further information on the survey fields. 

3. PHOTOMETRIC REDSHIFTS 

Photometric redshifts were derived using the method pre- 
sented in Bender et al. (2001). This method is a template match- 
ing algorithm rooted in Bayesian statistics. It closely resembles 
the method presented by Bemtez (2000). The templates are de- 
rived by fitting stellar population models of Maraston (1998) to 
combined broad-band energy distributions of MUNICS galax- 
ies with spectroscopic redshifts. In this way, representative 
galaxy templates of mixed stellar populations (variable age, 
metallicity, and dust extinction) optimized for the MUNICS 
dataset are obtained. 

The use of photometric redshifts can introduce systematic er- 
rors in the derived galaxy distances. To investigate these errors, 
we performed the following tests: 1) Templates were derived 
for a subsample of objects with spectroscopic redshifts and their 
suitability for the whole galaxy sample was verified using the 
remaining objects with spectroscopy. 2) Using subsets of the 
template library, we estimated the errors in the luminosity and 



mass functions introduced by a finite (and incomplete) set of 
templates. These were included in our error budget in the anal- 
ysis. 3) Monte-Carlo simulations were performed to investigate 
the influence of photometric redshift errors on the luminosity 
function (LF) and the galaxy number densities. The shape of 
the LF is not significantly biased as long as the photometric red- 
shifts scatter symmetrically around the true redshifts (see also 
Subbarao et al. 1996) and their errors are smaller than the bin 
size in z over which the luminosity function is averaged. 

Fig. 1 shows the comparison of photometric and spectro- 
scopic redshifts for 310 objects wifliin 5 MUNICS fields. The 
typical scatter in the relative redshift error Az/(1 +z) is 0.075. 
The mean bias is negligible. The distribution of photometric 
redshifts peaks around z ~ 0.5. and has a tail extending to z « 3. 

To further strengthen the confidence in the photometric red- 
shifts we show the restframe ^-band luminosity function (LF) 
obtained from these data in Fig. 2. 

The LF shows no significant evolution up to a redshift of 1 .2, 
in agreement with the finding of Cowie et al. (1996) from a 
much smaller spectroscopic sample. 

4. THE INTEGRATED STELLAR MASS FUNCTION 

The integrated stellar mass function > ^n^), the co- 
moving number density of objects having stellar mass exceed- 
ing is computed using the Vmax formalism to account for 
the fact that some fainter galaxies are not visible in the whole 
survey volume. Each galaxy in a given redshift bin [z;,Z;,) con- 
tributes to the number density an amount inversely proportional 
to the volume in which the galaxy is detectable in the survey: 

V, = da —dz, (1) 

Jz, dz 

where dV /dz is the comoving volume element, dQ. is the sur- 
vey area, Znua is the maximum redshift at which galaxy / having 
absolute magnitude ^ is still detectable given the limiting 
apparent magnitude of the survey and the galaxy's SED (the 
best-fit SED from the photometric redshift determination). 

Additionally, the contribution of each galaxy i is weighted by 
the inverse of the detection probabihty, P{m^ ■), where we as- 
sume that the detection probability is independent of the galaxy 
type and can be approximated by that of point-like sources. We 
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Fig. 2. — The restframe /T-band luminosity function derived from the MUNICS data in four redshift bins spanning 0.4 < z < 1.2. The dotted and the dashed curve 
are the z = LF by Loveday (2000) and Kochanek et al. (2001), respectively. 



use only objects with P(ot^ ,) > 0.75, such that this correction 
is always small. We have checked that this correction does not 
bias our results by comparing to what we get for higher com- 
pleteness Umits. 

To compute the stellar mass of a galaxy, we use an approach 
which maximizes the stellar mass for any /T-band luminosity 
at any redshift. Noting that ^ jL^ is a monotonically rising 
function of age for Simple Stellar Populations (SSPs), we find 
that the Ukely upper limit for jL^^ is the mass-to-Ught ratio 
of a SSP which is as old as the universe at the galaxy's redshift. 
This is the most extreme case of passive luminosity evolution 
(PLE) one can adopt. It corresponds to a situation where all 
massive galaxies would be of either elliptical, SO, or Sa type. 
We take the mass-to-light ratios from the SSP models published 
by Maraston (1998), using a Salpeter IMF. Similar dependen- 
cies on age are obtained from the models of Worthey (1994) 
and Bruzual & Chariot (1993) although the absolute values of 
M jL^ vary somewhat. The M jL^ which we obtain with our 
cosmological parameters at z ~ 0.5 are approximately consis- 
tent with ^ jL^ of local galaxies (Bell & de Jong 2001). To 
obtain a more reahstic estimate of ^ jL^, we used our VRIJK 
color information and the photometric redshift to fit the age and 
SFR of each galaxy using a grid of composite stellar popula- 
tions with exponential star formation timescales ranging from 
1 to 8 Gyr We obtain the same slope in the evolution of the av- 
erage jh^ with redshift as the one obtained under the PLE 
assumption. This modelhng will be described in detail in a 
forthcoming paper 

The mass in solar units is given by 

jt=[^\Q-^mK'M^,K)^ (2) 

using M,; ^ = 3.33. 

The number density in each redshift bin, ni^M > ^^^,z), is 
finally computed by summing over all galaxies in the bin whose 
stellar mass is exceeding ^]^^, 

n{^>J^,^,z)=Z^^y (3) 

The resulting integrated mass functions for ^^aa = 2 x 
lO'O/i-^^Q, ^li^ = 5 X lO'O/z-^^Q, and ^^^^ = 1 x 
IQ^^h^^^Q are shown in Fig. 3 along with the integrated lu- 
minosity functions for comparison. The mean values of^/Lj. 



in the maximum PLE model in the four redshift bins are 0.99, 

0.83, 0.73, and 0.65, as computed from the look-back time in 
our cosmology. With these mean values the mass limits cor- 
respond to absolute /T-band magnitudes of —22.43, —22.63, 
-22.77, and -22.90, respectively, for = 2 x \Q^^h-^J(.y 
For ^lij^ = 5 x \0^^h-^J(Q the numbers are -23.42, -23.62, 
-23.76, and -23.89. Finally, for J^^-^ = 1 x lO'^/z-^^g we 
have -24.18, -24.38, -24.51, and -24.64 (magnitudes with 
respect to h=\). 

5. DISCUSSION 

Fig. 3 compares the evolution of the integrated luminosity to 
the integrated mass. It is evident that the number density of lu- 
minous /T-band selected galaxies does not evolve significantly 
(given our uncertainties) to z = 1.2. However, because of the 
inevitable evolution of the mass-to-Ught ratio with z, the num- 
ber density of massive systems does change. Transforming lu- 
minosities into masses with our maximum PLE scheme yields 
a roughly constant number density for our lowest mass limit, 

2 x lO^^h^^^Q, and a decrease of the number density with red- 
shift by a factor of ~ 3 for a mass Umitof 5 X and 
by a factor ~ 6 for objects more massive than 1 x IQ^^h^^^Q. 
As the true jL^ at high redshift will most likely be lower 
than in our maximum PLE model, the true number densities are 
likely to decrease more rapidly with redshift. 

The steepening of the curves with increasing limiting mass 
in the maximum PLE case (despite them all having the same 
mass-to-Ught ratios at any given redshift) is due to the invari- 
ance of the LF with redshift and its steepness at the bright end. 
At increasing limiting mass, one is moving down the steepen- 
ing bright end of the LF, so that the same change in the mass- 
to-light ratio yields a higher change in the number density. 

The main uncertainty in this conclusion is still the field to 
field variation, in spite of the relatively large area surveyed, fol- 
lowed by the choice of SED templates used in the photometric 
redshift code (see above). The size of the open symbols in Fig. 

3 represents our estimate of the total uncertainty of the mean 
values. If we assume a Gould IMF (Gould et al. 1998) instead 
of a Salpeter IMF, the evolving ^/L^ curve becomes lower 
in its normalization as the mass-to-light ratio becomes smaller 
due to the reduced number of low-mass stars. The slope does 
not change significantly. 
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Fig. 3. — Comoving number density of objects having rest frame ^-band luminosities exceeding —22.42 + 51og/i (2 x 10'"/i L-.), — 23.42 + 51ogA 
(2 X lO'^/i^^Lg), and -24.17 + 51ogft (2 x lO'^h-^L,) (upper panels) and comoving number density of objects having stellar masses exceeding ^y^^ = 
2 X lO'^/j-^^g, = 5 X lO'^A^^^Q, and = 1 x 10"ft-2^Q (integrated stellar mass functions; lower panels). Mass to light ratios are assigned to 

maximize the stellar mass at a given luminosity (see text), and thus are likely upper limits. The solid points denote the values measured separately in each survey 
field, the open circles denote the mean values over the whole survey area. The size of the open circles is chosen to represent our estimate of the total uncertainty in 
the mean values. 



The observed density evolution as a function of mass is 
qualitatively consistent with the expectation from hierarchical 
galaxy formation models. Most rapid evolution is predicted for 
the number density of the most massive galaxies while the num- 
ber density of L*-galaxies should evolve much less. E.g. Baugh 
et al. (1998) predict that the number density of galaxies of a 
stellar mass of IO^^/z^^^q decreases by a factor of ~ 3.1 over 
redshift range 0.4 < z < 1 .2 (for the cosmological parameters as 
used here). Though this agreement is encouraging, both more 
elaborated models and improved sets of data are required. The 
latter can be obtained by larger and deeper samples, and more 
reaUstic estimates of jL^. 



We would Uke to thank the Calar Alto staff for their long- 
standing support during many observing runs over the last 
five years. This work was partly supported by the Deutsche 
Forschungsgemeinschaft, grant SFB 375 "Astroteilchenphysik" 
and the German Federal Ministry of Education and Research 
(BMBF), grant 05 AV9WM1/2. ND would like to thank S. 
White and G. Kauffmann for helpful discussions. GJH ac- 
knowledges support by the Texas Advanced Research Program 
under Grant No. 009658-0710-1999. CMdO acknowledges 
support from FAPESP (Fundagao de Amparo a Pesquisa do Es- 
tado de Sao Paulo) and the Alexander von Humboldt Founda- 
tion. 



REFERENCES 



Baugh, C. M., Cole, S., Frenk, C. S., & Lacey, C. G. 1998, ApJ, 498, 504 
Bell, E. R, & de Jong, R. S. 2001, ApJ, 550, 212 

Bender, R., et al. 2001, in ESO/ECF/STScI Workshop on Deep Fields, ed. 

S. Christiani (Berlin: Springer), 327 
Bender, R., Ziegler, B., & Bruzual, G. 1996, ApJ, 463, L51 
Bem'tez, N. 2000, ApJ, 536, 571 
Brinchmann, J., & Ellis, R. S. 2000, ApJ, 536, L77 
Bruzual, G. A., & Chariot, S. 1993, ApJ, 405, 538 

Cole, S., Lacey, C. G., Baugh, C. M., & Frenk, C. S. 2000, MNRAS, 319, 168 
Cowie, L. L., Songaila, A., Hu, E. M., & Cohen, J. G. 1996, AJ, 112, 839 
Drory, N., Feulner, G., Bender, R., Botzler, C. S., Hopp, U., Maraston, C, 

Mendes de Oliveira, C, & Snigula, J. 2001, MNRAS, 325, 550 
Ellis, R. S., Colless, M., Broadhurst, T, Heyl, J., & Glazebrook, K. 1996, 

MNRAS, 280, 235 
Gould, A., Flynn, C, & Bahcall, J. N. 1998, ApJ, 503, 798 
Kauffmann, G., & Chariot, S. 1998, MNRAS, 297, L23 
Kauffmann, G., White, S. D. M., & Guiderdoni, B. 1993, MNRAS, 264, 201 
Kelson, D. D., Illingworth, G. D., Franx, M., & van Dokkum, P. G. 2001, ApJ, 

552, L17 



Kochanek, C. S., et al. 2001, ApJ, submitted 

Lilly, S. J., Le Fevre, O., Crampton, D., Hammer, R, & Tresse, L. 1995, ApJ, 

455, 50 

Lin, H., Yee, H. K. C, Carlberg, R. G., & Ellingson, E. 1997, ApJ, 475, 494 
Lin, H., Yee, H. K. C, Carlberg, R. G., Morris, S. L., Sawicki, M., Patton, 

D. R., Wirth, G., & Shepherd, C. W. 1999, ApJ, 518, 533 
Loveday, J. 2000, MNRAS, 312, 557 
Maraston, C. 1998, MNRAS, 300, 872 
Rix, H., & Rieke, M. J. 1993, ApJ, 418, 123 
Somerville, R. S., & Primack, J. R. 1999, MNRAS, 310, 1087 
Subbarao, M. U., Connolly, A. J., Szalay, A. S., & Koo, D. C. 1996, AJ, 112, 

929 

van Dokkum, P. G., Franx, M., Kelson, D. D., & Illingworth, G. D. 2001, ApJ, 
553, L39 

Vogt, N. P, Forbes, D. A., Phillips, A. C, Gronwall, C, Faber, S. M., 

Illingworth, G. D., & Koo, D. C. 1996, ApJ, 465, L15 
Worthey, G. 1994, ApJS, 95, 107 



